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Abstract. We present the first results from the Complete Optical and Radio Absorption Line System (CORALS) 
survey. We have compiled a homogeneous sample of radio-selected QSOs from the Parkes Catalogue and searched 
for damped Lyman alpha systems (DLAs) towards every target, irrespective of its optical magnitude. This ap- 
proach circumvents selection effects - particularly from intervening dust - which have long been suspected to 
affect DLA surveys in optically-selected, magnitude-limited QSO samples. The CORALS data set consists of 66 
z cm > 2.2 QSOs in which 22 DLAs with absorption redshifts 1.8 < z a bs < z C m have been identified over a total 
redshift interval = 55.46 . Three of the DLAs are classified as 'associated' systems with z a bs ~ -z cm ; of the 19 
intervening DLAs, 17 are new discoveries. In this first paper of the CORALS series we describe the sample, present 
intermediate resolution spectroscopy and determine the population statistics of DLAs. We deduce a value of the 
neutral gas mass density traced by DLAs (expressed as a fraction of the closure density) log QohAh = — 2.591q'24, 
and a number density of DLAs per unit redshift n{z) = 0.3ltg;g§, both at a mean redshift (z) — 2.37. Both values 
are only marginally higher than those measured in optically selected samples of QSOs. Taking into account the er- 
rors, we conclude that dust-induced bias in previous surveys may have led to an underestimate of these quantities 
by at most a factor of two. While n(z) is greater in fainter (B > 20) QSOs, the effect is only at the ~ la level and 
we have not uncovered a previously unrecognised population of high column density (iV(H I)> 10 21 cm" 2 ) DLAs 
in front of faint QSOs. These conclusions are tentative because of the limited size of our data set; in particular the 
distribution of column densities is poorly sampled at the high end where a much larger survey of radio-selected 
QSOs is required the improve the statistics. 
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1. Introduction 

The presence of intervening gas in the sightline to a high 
redshift QSO will cause absorption lines to be superim- 

Send offprint requests to: S. Ellison 

* The work presented here is based in part on data ob- 
tained with the ESO facilities on La Silla (EFOSC/3.6-m) and 
Paranal (FORS1/UT1). 



posed on the quasar's continuum. This relatively simple 
physical process allows us to study in detail the properties 
of gas clouds in front of QSOs, whether in the intergalac- 
tic medium or in the galaxies giving rise to high column 
density damped Lyman alpha Systems (DLAs). DLAs are 
one of the key components of the universe at high redshift 
because, although relatively rare, they account for most 
of the neutral gas available for star formation. 
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Measuring the redshift evolution of 12dla (the total 
amount of neutral gas traced by DLAs expressed as a 
fraction of the closure density) has traditionally been seen 
as a tool for probing the history of assembling galaxies 
and measuring the rate at which they convert gas into 
stars (e.g. Lanzetta, Wolfe, & Turnshek 1995; Pei & Fall 
1995; Storrie-Lombardi, McMahon & Irwin 1996; Storrie- 
Lombardi & Wolfe 2000; Rao & Turnshek 2000). Similarly, 
chemical abundances in DLAs have been used to trace the 
metallicity evolution of galaxies over most of the Hubble 
time (e.g. Pettini et al. 1997; Pettini et al. 1999; Prochaska 
& Wolfe 1999, 2000). However, the interpretation of both 
DLA abundances and S2dla as indicators of galaxy evo- 
lution is pivotal upon the assumption that DLAs are rep- 
resentative of the bulk of 'normal' galaxies at each red- 
shift sampled. Therefore, it is rather surprising that, at 
least for current samples of DLAs, no strong evolution of 
either SIdla or metallicity (Zdla) has been revealed by 
the above studies from z ~ 3.5 down to z ~ 0.5 . One 
possible explanation is that DLAs are representative of 
only a particular evolutionary stage in a galaxy's lifetime. 
Alternatively, current DLA samples may be incomplete 
due to a selection bias In this paper, we specifically inves- 
tigate the possibility that dust extinction may have pro- 
duced incomplete samples of DLAs in previous surveys 
based on optically-selected QSOs. 

DLAs appear to remain relatively metal-poor (Zdla ~ 

— jq-Z©) from z ~ 3.5 to z ~ 0.5 (Pettini et al. 1999; 
Prochaska & Wolfe 2000). However, it has been noted by 
several authors that the absorbers with the highest val- 
ues of hydrogen column density are all of low metallic- 
ity and dust content (e.g. Pettini et al. 1997; Boisse et 
al. 1998; Pettini et al. 1999; Prantzos and Boissier 2000; 
Savaglio 2001). Whilst relatively metal-rich damped sys- 
tems do exist, they tend to have values of N(H I) to- 
wards the lower end of the range considered for DLAs, 
N(R I)> 2 x 10 20 cm -2 . Boisse et al. (1998) and Prantzos 
& Boissier (2000) in particular have pointed out an appar- 
ent anti-correlation between 7V(H I) and metallicity. Since 
the census of metal abundances (using the column-density 
weighted average) is dominated by the highest A^(H I) ab- 
sorbers, such an anti-correlation could explain the lack of 
metallicity evolution in the present samples of DLAs. This 
anti-correlation between iV(H I) and metallicity may re- 
flect a bias against high column density, metal-rich DLAs 
with sufficient dust to obscure background QSOs from the 
eyes of magnitude-limited optical surveys. 

Further evidence that optical magnitude-limited QSO 
samples may be incomplete due to dust bias comes from 
consideration of the spectral energy distribution of quasars 
with and without DLAs. Fall, Pei, & McMahon (1989) and 
Pei, Fall, & Bechtold (1991) found evidence that QSOs 
with DLAs have systematically steeper continuum slopes 
than quasars with no damped absorber. A similar conclu- 
sion was reached by Carilli et al. (1998) who found that 
a high fraction of their 'red' QSO sample had associated 
21cm absorbers compared with QSOs selected solely for 
the presence of Mg II absorbers. By modelling the extinc- 



tion effect of dust, Pei & Fall (1995) estimated that at 
z = 3 between 27% and 44% of QSOs and 23%-38% of 
DLAs may be missing from i3-band selected quasar sam- 
ples. 

The objective of the Complete Optical and Radio 
Absorption Line System (CORALS) survey presented here 
is simple and straightforward; compile a homogeneous 
sample of radio-selected QSOs and obtain medium resolu- 
tion spectra of every target, regardless of optical magni- 
tude. In this way it should be possible to determine quan- 
titatively the severity of the dust bias implied by the ob- 
servations reviewed above. The present paper is the first 
in the CORALS series. Here we define the sample (§U), 
describe the observations and present the QSO spectra 
(§||), identify the DLAs (§||), and calculate their number 
density and 12dla (®- 

We adopt an f^f = 1.0,12a = cosmology and use 
Ho = 65 km s _1 Mpc _1 throughout the paper. 

2. Previous DLA Surveys 

It is useful at this point to review some of the previous 
surveys for damped Lyman a systems and summarise their 
sample definitions and major findings. Relevant data are 
collected in Table |l|. 

The first major survey for DLAs was conducted by 
Wolfe et al. (1986) (WTSC86) from the relatively bright 
Lick QSO sample. The strategy adopted by those authors 
was to obtain low resolution spectra for a large number 
of QSOs in order to identify absorption systems above a 
certain equivalent width limit (in this case Wo > 5 A). 
These candidate DLAs were then followed-up with 2 A 
resolution spectroscopy revealing that approximately 50% 
were indeed damped, the remainder being blends of lower 
column density lines. This strategy of pre-selecting DLA 
candidates from low resolution spectra, based on rest- 
frame equivalent width, has come to characterise surveys 
in subsequent years. The approach allows a large sample of 
QSOs to be considered, maximising efficient use of observ- 
ing time. The WTSC86 survey represents a landmark in 
absorption line studies, being the first major compilation 
of spectra and line lists for DLAs. One of the main results 
was that the incidence of DLAs was found to be greater 
than expected if the absorption were due to galaxies with 
the same cross section for H I absorption as present-day 
spirals (Wolfe 1988). In addition, the procedures for iden- 
tifying DLAs were established and the N(H I) definition 
of a damped system set by this work; N(H I)> 2 x 10 20 
cm~ 2 has been upheld in all subsequent studies. 

A survey of similar size to the WTSC86 sample 
was conducted a few years later by Sargent, Steidel k. 
Boksenberg (1989) (SSB89), although the initial focus of 
this work was on Lyman limit systems (LLSs). All of the 
quasars in the SSB89 sample were observed at a resolu- 
tion of 4 A, considerably higher than the first pass made 
by Wolfe et al. (1986). However, SSB89 estimate that their 
DLA sample is probably only complete for systems with 
N(H I)> 10 21 cm -2 , although several lower column den- 
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Survey reference Abbreviation No. QSOs Supplemented Limiting Mag. Telescope Resolution 

with (Band) 



Wolfe et al. (1986) 


WTSC86 


68 




18.5(10 


Lick 


io A 


Boksenberg (1989) 


GQT3CQ 


oo 




N/A a 

IN / /\ 


jnaiti u-m 


4 A 


ijaiizeLia eii txL. ^lyyij 


T ANQ1 


^7 
o t 


oo-Doy 


Ijliei dLUI c 

search 


I— 1 did til || t i | 1 o tl 

jnaie o-m diiu ±ja» 
Campanas 2.5-m 


A - (\ k 


Lanzetta, Wolfe & 
Turnshek (1995) 


LWT95 


260 


LAN91 


20(V) b 


IUE 


7- 10 A 


Wolfe et al. (1995) 


WLFC95 


228 


LWT95 


18.75(B) 


MMT and Las 
Campanas 2.5-m 


6 - 10 A 


Storrie-Lombardi 
& Wolfe (2000) 


SLW00 


40 


Many incl. WTSC86, 
LWT95, WLFC95 


19.5(i?) 6 


Keck, A AT 
Lick 


2 A 


Peroux et al. (2001a) 


PER01 


66 


SLW00 


20.5(R) b 


4-m CTIO 
WHT 


5 A 



Table 1. Summary of past surveys and their sample definitions. Column 3 tabulates the number of QSOs new to each 
survey Therefore, if a particular survey has been supplemented with a previous sample, the total number of QSOs 
is the sum of the individual entries in column 3. Notes: a - No broad-band limit is quoted for SSB89, only a derived 
apparent magnitude from the continuum flux at Ao = 1450 A. b - This is the faintest magnitude in this sample, but 
it does not indicate a completeness limit. 



sity lines were identified. Combining the SSB89 dataset 
with a similar number of QSOs found in the literature, 
Lanzetta et al. (1991) adopted the strategy of WTSC86 
in making a 5 A equivalent width cut to select 89 DLA 
candidates from their sample of 101 QSOs. 

From follow-up spectroscopy at 2 - 3 A resolution of 
these candidates, Lanzetta et al. (1991) produced the first 
thorough statistical analysis of damped system number 
density, the mass density of neutral gas in DLAs and 
their clustering properties. A few years later, Lanzetta et 
al. (1995) extended their work to include DLAs at lower 
rcdshift by exploiting ultraviolet data obtained from the 
International Ultraviolet Explorer (IUE) satellite to gain 
one of the first glimpses into the z a ^ s < 1.6 absorber pop- 
ulation. These authors found that I^dla decreases signif- 
icantly from z ~ 3.5 to z ~ 0.01, and interpreted this 
rcdshift evolution as evidence for the consumption of H I 
gas by star formation. 

Wolfe et al. (1995) (WLFC95) used the Large 
Bright Quasar Survey (LBQS, Hewett, Foltz & Chaffee 
1995 and references therein) to search for DLAs with 
1.6 £ z abs £ 3.0. A total of 59 DLA candidates out of 228 
spectra were pre-selected as having Wo > 5 A. At the 
time of publishing their paper, Wolfe et al. (1995) had con- 
firmed the identification of 13 DLAs with N(H i)> 2x 10 20 
cm -2 f rom 15 candidates with Wq > 9 A. In addition, 
there were 8 DLA candidates whose equivalent widths ex- 
ceeded 10 A and were therefore considered highly likely to 
be damped systems. All of these 8 candidates have subse- 



quently been confirmed as DLAs with intermediate resolu- 
tion spectroscopy (Storrie-Lombardi & Wolfe 2000). In ad- 
dition to the LBQS sample, WLFC95 constructed a 'sta- 
tistical sample' from the literature consisting of 80 DLAs. 
From this they confirmed the coincidence which had pre- 
viously been noted (e.g. Lanzetta 1993) that the mass of 
H I in DLAs at z ~ 3 is similar to the total luminous 
mass in stars today. Expressing both as fractions of the 
closure density, Q D la(^ = 3.25) = (5.1 ± 1.7) x lO" 3 /^ 1 
and fl*(z = 0) ~ 5 x lO" 3 ^ 1 (e.g. Pagel 2000). 

In order to extend studies of DLAs to lower red- 
shifts, Rao & Turnshek (2000) have recently published re- 
sults from an observing campaign with the Hubble Space 
Telescope (HST). Space-based telescopes are required in 
order to detect low z DLAs since the Lyman a signa- 
ture at z a b s £1.6 falls below the atmospheric cut-off. Since 
such satellite resources are limited, a pre-selection based 
on Mg II and Fc II line strengths determined from ground- 
based observations has permitted an efficient screening of 
the initial QSO sample. A total of 12 DLAs with z a b s < 1.5 
was found, thereby significantly improving the statistics 
in this redshift regime. The puzzling result obtained from 
these new data is that, contrary to previous indications, 
^dla appears to remain approximately constant over all 
redshifts sampled, from z = 3.5 to 0.1 . Using HST archival 
spectra, Churchill (2001) has also argued that the num- 
ber density of DLAs implied from the incidence of Mg II 
systems remains unaltered down to z a b s ~ 0.05. 
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Two recent surveys have pushed the bounds to higher 
redshifts. Storrie-Lombardi & Wolfe (2000) and Peroux 
et al. (2001a) have recently identified 19 and 15 DLAs 
respectively at z a b s > 3.5, including the highest redshift 
damped absorber known to date (z a bs = 4.466, Peroux et 
al. 2001a; Dessauges-Zavadsky et al. 2001). Again, there 
is only marginal evidence for evolution; between z = 3.5 
and 4.5 I^dla is lower than at more recent epochs, but the 
effect is only at the ~ 2a level. 

The intention of this (non-exhaustive) review of previ- 
ous DLA surveys is to illustrate the difficulty of extract- 
ing unbiased statistics from samples based on inhomoge- 
neous datasets. Many of these surveys have built on exist- 
ing samples, taken additional targets from the literature 
and sometimes have insufficient resolution to reliably de- 
termine whether an absorber is damped. The advantage 
of this approach has been largely a statistical one. The ag- 
gregate samples often contain several hundred QSOs and 
such a large number is indeed required in order to deter- 
mine meaningful statistics because of the rarity of damped 
systems. The purpose of CORALS is to provide the first 
complete and homogeneous survey for DLAs free from any 
bias that may be inherent in magnitude-limited QSO sam- 
ples. 



Telescope Date Resolution No. of 

(no. nights) (at 4000 A) QSOs observed 



ESO 3.6-m 


Sept. 1998 (2) 


6 - 8 A 


16 


ESO 3.6-m 


Feb. 1999 (2) 


7 A 


12 


AAT 


Dec. 1998 (3 a ) 


3 A 


11 


AAT 


Apr. 1999 (2.5 b ) 


3 A 


10 


AAT 


Oct. 1999 (2) 


3 A 


12 


VLT 


Oct. 2000 (0.5 C ) 


4.5 A 


9 


VLT 


Mar. 2001 (2) 


4.5 A 


10 


VLT 


June 2001 (0.5 d ) 


4.5 A 


1 



Table 2. Summary of observations CORALS survey. 
Many of the bright targets were observed first with low 
resolution at the ESO 3.6-m and QSOs with strong ab- 
sorbers followed-up with AAT observations, whilst others 
were observed directly with the AAT. All faint targets 
(B > 20) were observed only with the VLT. Notes: a - 
One night lost to bad weather. b - The April AAT run 
consisted of 5 half nights. c - 6 hours of service time. d - 
6 hours of Director's discretionary time. 



3. The CORALS Survey Sample 

We have based our new survey for DLAs on the radio- 
selected Parkes quarter- Jansky sample (Shaver et al. 1996; 
Jackson et al. in prep.; Hook et al. in prep.). This com- 
plete sample of flat-spectrum radio sources was taken from 
the Parkes Catalogue (Wright & Otrupcek 1990), which 
consists of radio and optical data for almost 8300 radio 
sources for the sky south of declination +27deg. Over 
most of its area, it is complete to S2.7 GHz = 0.25 Jy. 
The Parkes quarter- Jansky sample consists of all (878 in 
total) flat-spectrum (a > —0.4, measured at 2.7 and 5.0 
GHz) sources with declinations between +2.5° and —80°, 
excluding low galactic latitudes (| b |< 10°) and regions 
around the Magellanic Clouds. The sample is complete to 
5*2.7 GHz = 0.25 Jy over 4.39 sr and to 0.60 Jy over a 
further 1.16 sr. 

Using this Parkes compilation of radio sources as the 
parent sample, optical identification and i?-band magni- 
tudes for all targets with S > —45° were achieved using the 
COSMOS Southern Sky Catalogue (Drinkwater, Barnes 
& Ellison 1995) in the first instance or imaging at the 
ESO 3.6-m telescope at La Silla for the faintest QSOs. 
Therefore, optical counterparts were identified for every 
radio source in the Parkes compilation. From this complete 
sample of 878 sources with no optical magnitude limit or 
other selection bias, low resolution spectra (FWHM = 12 
- 14 A) were obtained for the 442 stellar identifications 
(QSOs and BL Lacs) with the EFOSC on the ESO 3.6-m 
to determine redshifts. We note that since compiling the 
CORALS sample, the catalogue of PKS sources has been 
revised slightly to include some extra targets (Jackson et 
al. in prep.; Hook et al. in prep.). However, it is impor- 



tant to stress that this will not affect the results presented 
here, since the important factor in our sample definition 
is that it is optically complete. 

In our ground-based search for DLAs in this sample 
we restricted ourselves to QSOs with z cm > 2.2, so as 
to be able to record a sufficient portion of the spectrum 
blueward of the Lya emission down to the onset of the 
atmospheric cut-off. This left us with a final sample of 66 
QSOs in which we could search for DLAs in the range 
1.8 < z a bs < z cm with sufficient sensitivity to measure 
reliable values of the column density N(R I).Q 

4. Observations and Spectral Analysis 

Before commencing our observing campaign, we first 
searched the literature to determine which QSOs had al- 
ready been observed at sufficiently high spectral resolu- 
tion, S/N and over the correct wavelength range to match 
the rest of our DLA survey. The eight QSOs for which 
adequate spectra were found in the literature were not 
re-observed. The remaining 58 QSOs were subject to an 
extensive observing campaign that spanned five semesters 
on three different telescopes; relevant details are presented 
in Table ||. References for the eight QSO spectra observed 
by others are given in Table |. 

Our observing strategy has been to obtain 'snapshot' 
spectra of as many bright (typically B < 20) QSOs as 
possible at low resolution with the ESO 3.6-m in order to 
pre-select targets with candidate DLAs for higher resolu- 
tion follow-up at the AAT. Targets fainter than this limit 

1 One of the observed QSOs actually has z em < 2.2 due to 
an incorrect redshift in our original list, see Table 0. 
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QSO 




B Mag 


DLA? 


N(B I) (10 20 cm- 2 ) 




Re 


B0017-307 


2.666 


19.0 


N 








B0039-407 


2.478 


18.5 


N 








B0104-275 


2.492 


18.5 


N 








B0113-283 


2.555 


19.0 


N 








B0122-005 


2.280 


18.5 


N 








B0244-128 


2.201 


18.5 


N 








B0256-393 


3.449 


19.6 


N 








B0325-222 


2.220 


19.0 


N 








B0329-255 


2.685 


17.1 


N 






1 


B0335-122 


3.442 


21.5 


Y 


6.0 


3.178 




B0347-211 


2.944 


21.1 


Y 


2.0 


1.947 




B0405-331 


2.570 


19.0 


Y 


4.0 


2.570" 




B0420+022 


2.277 


19.5 


N 








B0422-389 


2.346 


18.0 


N 








B0432-440 


2.649 


19.6 


Y 


6.0 


2.297 




B0434-188 


2.702 


20.0 


N 








B0438-436 


2.863 


19.5 


Y 


6.0 


2.347 




B0451-282 


2.560 


19.0 


N 








B0458-020 


2.286 


20.0 


Y 


45.0 


2.039 


2 


B0528-250 


2.765 


19.0 


Y 


5.6 


2.141 


3 








Y 


15.8 


2.811" 


4 


B0537-286 


3.110 


20.0 


Y 


2.0 


2.974 




B0601-172 


2.711 


20.0 


N 








B0610-436 


3.461 


19.0 


N 






5 


B0819-032 


2.352 


18.2 


N 








B0834-201 


2.752 


19.0 


N 








B0913+003 


3.074 


21.7 


Y 


5.5 


2.744 




B0919-260 


2.300 


19.0 


N 








B0933-333 


2.906 


20.0 


Y 


3.0 


2.682 




B1010-427 


2.954 


17.5 


N 








B1055-301 


2.523 


19.5 


Y 


35.0 


1.904 




B1136-156 


2.625 


20.0 


N 








B1147-192 


2.489 


19.4 


N 








B1149-084 


2.370 


18.5 


N 








B1228-113 


3.528 


22.0 


Y 


4.0 


2.193 




B1228-310 


2.276 


19.0 


N 








B1230-101 


2.394 


19.8 


Y 


3.0 


1.931 




B1251-407 


4.464 


23.7 


Y 


4.0 


3.533 










Y 


2.0 


3.752 




B1256-243 


2.263 


19.5 


N 








B1318-263 6 


2.027 


20.4 


N 








B1351-018 


3.710 


20.9 


N 








B1354~107 


3.006 


19.2 


Y 


2.5 


2.501 










Y 


6.0 


2.966 a 




B1402-012 


2.518 


18.2 


N 








B1406-267 


2.430 


21.8 


N 








B1418-064 


3.689 


18.5 


Y 


2.5 


3.449 




B1430-178 


2.331 


19.0 


N 








B1535+004 


3.497 


24.1 


N 








B1556-245 


2.813 


18.5 


N 








B1635-035 


2.871 


21.8 


N 









Table 3. Complete list of all CORALS QSOs and the tally of DLAs identified. References refer to QSOs for which 
optical spectroscopy had already been obtained prior to our survey, as follows: 1 - Lanzetta et al. (1991); 2 - Pettini 
et al. (1997); 3 - Ledoux et al. (1998); 4 - Lu et al. (1996); 5 - M. Giavalisco, private communication; 6 - Tytler, Fan 
and Buries (1996); 7 - Songaila (1998); 8 - Turnshck et al. (1989) Notes: a - z abs ~ z cm ; a - The original rcdshift of 
this target was incorrect in our list, the VLT spectrum confirms the true redshift which is less than the survey cut-off 

&f 2em — 2.2. 
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QSO z em B Mag DLA? N(R I) (10 20 cm -2 ) z abs Ref 



B17014 


-016 


2.842 


21.7 


N 


B1705-f 


-018 


2.575 


18.9 


N 


B1937- 


-101 


3.780 


19.0 


N 


B2000- 


-330 


3.780 


18.5 


N 


B2126- 


-158 


3.275 


17.5 


N 


B2149- 


-307 


2.330 


17.5 


N 


B2212- 


-299 


2.703 


17.8 


N 


B2215+020 


3.550 


21.5 


N 


B2224+006 


2.248 


21.7 


N 


B2245- 


-059 


3.295 


19.5 


N 


B2245- 


-328 


2.268 


16.5 


N 


B2256+017 


2.663 


19.0 


N 


B2311- 


-373 


2.476 


18.5 


Y 


B2314- 


-340 


3.100 


18.5 


N 


B2314- 


-409 


2.448 


18.0 


Y 










Y 


B2315- 


-172 


2.462 


19.5 


N 


B2325- 


-150 


2.465 


19.5 


N 


B2351- 


-154 


2.665 


19.0 


N 



3.0 

4.0 
2.0 



2.182 

1.857 
1.875 



Table 3. Continued 



were observed directly with the VLT. We briefly review 
the main characteristics of the observations. 

4.1. ESO 3.6-m Observations 

We obtained low resolution spectra of 28 QSOs using the 
EFOSC2 spectrograph on the ESO 3.6-m telescope during 
four nights between Sept. 1998 and Feb. 1999. The spec- 
tral resolution is approximately 7 A with a 1 arcsec slit 
and B600 grism, covering 3270 - 5240 A. For a few QSOs, 
we also took spectra using R600 grism with the similar 
spectral resolution to cover 4320 - 6360 A. The seeing 
conditions remained around 1.2 arcsec, although the sky 
was not photometric during all four nights. However, since 
we normalise our spectra to fit the DLA profile, the accu- 
rate flux calibration is not critical for the purpose of this 
paper. 

4.2. A AT Observations 

The RGO spectrograph was used with the TEK CCD, 
25cm camera and 600 V grating. The seeing at Siding 
Spring was highly variable through these observations, 
ranging from 0.8 to 2.0 arcsecs, often with large fluctu- 
ations during a given night. However, the slit width was 
fixed at 1.5 arcsec which projected onto 1.8 pixels to give 
a FWHM resolution of 3 A. The grating angle was chosen 
to cover the entire wavelength range of each QSO from at 
least 3400 A (corresponding to the wavelength of Lyman cv 
A1216 at z abs = 1.8) to z cm . The grating angle varied 
slightly for each run depending on the QSO subset to be 



observed, but was typically around 19.74°, corresponding 
to a central wavelength of A = 4000 A. At this setting, 
the wavelength range is ~ 3200 A - 4800 A. For the few 
bright, high redshift (z cm ^ 2.9) targets observed with the 
AAT, a second grating setting, with a central wavelength 
A = 5500 A, was required to cover the spectrum up to 
6300 A. 

4.3. VLT Observations 

The VLT observations were executed with FORS1 in a 
combination of service and visitor mode. The 600B grating 
was used for all targets, with additional wavelength cover- 
age provided by the 600R grating for the z om = 4.458 QSO 
B1251-407. A slit width of 0.7 arcsec provided a FWHM 
resolution of 4.5 A and a wavelength coverage of 3360 - 
5760 A. The 600R grating gave additional coverage (only 
required for the highest redshift QSO in our sample) over 
5200 - 7300 A with a resolution of 3.6 A. Despite the faint- 
ness of several of the VLT targets (down to B = 24) , all 
acquisitions could be executed in 'fast' mode and without 
blind offsets. 

4.4. Data Reduction and Column Density Fitting 

We applied the same reduction procedure to all of the 
data. The standard IRAF0 routines were used. First, all 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc, under coopera- 
tive agreement with the National Science Foundation. 



S. L. Ellison et al.: CORALS I: DLA number density and Odla- 



7 



of the images were trimmed and the bias level was sub- 
tracted using the over-scan regions. High S/N flat-field 
images were combined to a single image, which was then 
smoothed using a box median filter of 1 x 50 pixels. The 
orientation of the box median filter is chosen such that 
we preserve all of the variations along the dispersion axis. 
Then, the normalized 2-D flat-field image was used to re- 
move pixel-to-pixel variations in the quasar images. The 
task APALL was used to perform the optimal extraction 
of the 1-D spectra. Wavelength calibration images were 
typically taken both before and after each QSO image. 
The comparison lamps are CuAr at AAT, NcHc at ESO 
3.6m and NeHgCd on FORS1 at VLT. The dispersion so- 
lutions were obtained using a 4th order Legendre polyno- 
mial and the RMS errors of the fitting was less than 0.05A. 
Finally, we performed error-weighted summation of all of 
the wavelength calibrated spectra for each QSO. The final 
error array was the quadratic sum of the individual error 
spectra. 

The spectra of all 58 QSOs observed by us are pre- 
sented in Figure |]. Once extracted, the spectra were in- 
spected for the presence of DLAs - the Lyman a signature 
clearly visible as a broad, saturated absorption feature. If 
a DLA was identified, the spectrum was normalised by 
dividing through by the QSO continuum and then fitted 
with a Lyman a profile using the Starlink package D1PSO 
to determine the redshift and column density of the DLA. 
The complete list of CORALS QSOs and identified DLAs 
can be found in Table ||. 

5. Damped Lyman Alpha Systems in the 
CORALS Survey 

We adopt the usual definition of a DLA, i.e. N(H I) 
> 2 x 10 20 cm -2 , although our data are of sufficient res- 
olution and S/N to recognise and measure absorbers of 
somewhat lower 7V(H I) (Ellison 2000). From the observa- 
tions of the 66 CORALS QSOs, a total of 22 DLAs have 
been identified. Three of these have absorption redshifts 
similar to the emission redshift of the QSO. We follow the 
standard procedure of excluding DLAs within 3000 km s _1 
of the QSO redshift from our statistical analysis, in order 
to facilitate comparison with other surveys. However, we 
note that the z a b s ~ z CT[L DLAs are probably similar to in- 
tervening absorbers (M0ller, Warren & Fynbo 1998) and 
we defer the analysis of the 'associated' CORALS DLAs 
to a future paper (Ellison et al., in preparation). Of the 
19 intervening DLAs, two (B0458-020 and B0528-250a) 
were already known and have been extensively studied 
prior to our survey. Profile fits to all the others are shown 
in Figure ||. We now briefly discuss each DLA system. 

B0335— 122 Extended absorption around the DLA trough 
complicates the fitting of this system. The O I A1302 line 
is used to constrain redshift, but the error on the final fit 
is relatively high (6 ± 1 x 10 20 cm -2 ) due the critical, yet 
uncertain, continuum placement. 



B0347— 211 The combination of low absorption redshift 
and faint QSO magnitude have resulted in a low S/N spec- 
trum of this DLA and a relatively poor fit. The redshift 
is constrained by the Al II A 1670 line and the N(H I) 
determined by fitting a damped profile is in good agree- 
ment with the column density inferred from the equivalent 
width measurement. 

B0432— 440 Although the spectrum is noisy, the column 
density of this DLA is reasonably constrained by the base 
of its trough and the shape of its red wing. This system is 
well-fitted with a column density A^(H I) = 6 X 10 20 cm -2 . 

B0438-436 Both AAT and VLT spectra were obtained 
for this QSO. The best fit to the combined spectrum 
has A^(H I)= 6 x 10 20 cm -2 , although this is constrained 
mostly by the red wing since there is significant contam- 
inating absorption in the blue wing. The redshift of the 
Al II A1670 line (with rest-frame equivalent width Wq = 
0.63 A) agrees well with the redshift determined from the 
centre of the DLA trough. 

B0537— 286 A relatively high S/N spectrum and simple 
structure around the absorber permits a good fit to this 
DLA. 

B0913+003 An excellent fit to this DLA is facilitated by 
the clearly defined damping wings and lack of blending. 

B0933— 333 Despite blending with a weaker component 
blueward of the DLA, this system is reasonably fit with 
profile of iV(H I) = 3 x 10 20 cm" 2 . 

B1055— 301 This DLA has a very large equivalent width 
(Wo ~ 60 A) but is heavily blended with other absorption 
lines. The fit is only constrained by the base of the ab- 
sorption which is clearly saturated over 10 A in the rest 
frame. Higher order Lyman lines are not available to im- 
prove the decomposition of the H I cloud model, although 
several metal transitions are covered by the AAT spec- 
trum. Strong Si II A1526, Al II A1671 and Fe II A1608 are 
all observed with redshifts of 1.9037 and rest frame equiv- 
alent widths Wo = 1.2 A, 1.3 A and 0.8 A respectively. 
This provides strong support for the presence of a DLA 
at the position shown in Figure |[ However, the column 
density can only be constrained to within ±15%: A^(H I)= 
35 ± 5 x 10 20 cm" 2 . 

B1228— 113 Another low redshift system whose spec- 
trum has a low S/N, although the lack of strong nearby 
Lyman a forest lines results in an acceptable fit. 

B1230— 101 Constrained mostly by its fit to the red wing, 
this DLA at z a b s = 1.931 has several metal lines associated 
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with it. The AAT spectrum covers both Fe II A1608 (W 
= 520 mA) and Si II A1526 (W = 720 mA); the latter 
lies just blueward of a strong, resolved, C IV doublet at 
Zabs — 1.899. There is a second C IV system at A D b s ~ 4540 
A associated with the DLA itself. 

B1251— 407a, b As can be seen in Figure 0, this QSO has 
two prominent absorption lines at approximately 5510 and 
5780 A, corresponding to Lyman a at z a b s = 3.533 and 
3.752 respectively. The former is well reproduced by a 
damped profile with N(R I) = 4 x I0 20 cm~ 2 . The lat- 
ter has a large equivalent width (Wo = 13. 7 A), but also 
steep sides and cannot be fitted satisfactorily with a single 
absorption component. Closer inspection reveals structure 
within the core of the damped Lyman a line (see Figure 
2). We consider this feature to be a composite consisting 
of a DLA at z abs = 3.752 with N(R I) = 2 x 10 20 cm~ 2 
flanked by two lower column density components. This in- 
terpretation is supported by the presence of Si II A1526 
absorption at the same redshift as the DLA. In any case, 
both DLAs in this QSO are excluded from our discussion 
of the sample statistics below because we consider only the 
redshift interval 1.8 < z a b s < 3.5 (see §|J). The two other 
moderately large EW systems at A = 4640 and 4830 A 
respectively (see Figure 1), are not DLAs, but blends of 
lower column density lines. 

B1354— 107a Two strong absorption features are seen to- 
wards this QSO. Absorber 'a' has a redshift of z a b s = 2.501 
and is well fitted by a damped profile with N(H I) = 
2.5 x I0 20 cur 2 . Absorber 'b', with iV(H I)= 6 x I0 20 
cm~ 2 and z a b s = 2.966, is classified as a z a b s ~ z cm DLA 
and is not included in the present analysis. 



B1418-064 At z abs = 3.449, this is the highest red- 
shift DLA to be included in our statistical analysis; the 
Lyman a line is well reproduced by a damped profile with 
N(R I) = 2.5xl0 20 cm" 2 . Our limited spectral cover- 
age redward of Lyman a emission reveals associated Si II 
A1304 and O I A1302 absorption lines with W = 290 and 
390 mA respectively (although the O I transition is prob- 
ably blended with an unidentified line and the equivalent 
width measurement is therefore an overestimate). 



B2311-373 A column density of N(K I) = 3 xlO 20 cm" 2 
provides a satisfactory fit to the base and both wings of 
this DLA. Si II A1526 is the only metal transition identi- 
fied in the AAT spectrum with Wq = 500 mA. 

B2314-409a,b As is the case for B1055-301, there is ex- 
tended H I absorption in this spectrum. The spike at 
A ~ 1220 A in Fi gure 2 could be either noise or residual 
continuum flux between two closely spaced components. 
The latter interpretation is supported by the presence of 
two Si II A1526 absorption lines at redshifts z a b s = 1.857 



and 1.875 respectively. These redshifts match well the blue 
and red components of the Lyman a absorption feature, 
as indicated by the fit shown in Figure 2. The metal lines 
are stronger in the z a b s = 1.857 system, where we also 
detect Fe II A1608 and Al II A1671; the z abs = 1.875 com- 
ponent may well be a very low metallicity system given 
that Fe II A1608 has W < 90 mA. This is the second case 
of a multiple DLA, the first being a triple DLA (CTQ247) 
spread over ~ 6000 km s _1 , discovered by Lopez et al. 
(2000). 

6. DLA Statistics - Number Density and SIdla 

One of the key areas of interest in the study of DLAs is 
how the population evolves with time. If DLAs represent 
the bulk of the galaxy population at a given epoch, then 
studying their properties as a function of redshift will be 
a powerful method for tracing galaxy evolution. For ex- 
ample, there is strong evidence for redshift evolution in 
the column density distribution function, f(N), (Storrie- 
Lombardi & Wolfe 2000). However, since the number den- 
sity is well represented by a power law there appears to 
be no change in the product of space density and absorber 
cross section. As discussed in the Introduction, the evo- 
lution of f^DLA is still unclear, particularly at z a b s < 1.5 
where statistics are poor. 

Due to the relatively small size of the CORALS sur- 
vey, it is not possible to investigate the evolution of DLA 
statistics for this sample. However, since our main objec- 
tive is to ascertain whether or not a significant fraction of 
gas has gone undetected, it is sufficient for us to restrict 
our determination of Odla to the range 1.8 < z a b s < 3.5 
where there appears to be little evolution. This requires 
us to omit two DLAs from our sample (B1251— 407a and 
b), and restrict our statistical analysis to the remaining 
17 DLAs. 

6.1. DLA Number Density 

This is simply the total number of DLAs divided by the 
total intervening redshift interval covered, Az, given by 



(1) 



where the summation is over the n QSOs in the sample. 
For CORALS, z min = 1.80 and z max = 3.5 or z 30 oo (the 
redshift corresponding to a v — 3000 km s _1 ), whichever 
is smaller: 



^3000 




(b-l) 
(-6-1) 



- 1 



(2) 



3 We note that although most surveys calculate z m i n for each 
QSO on a case by case basis, depending on wavelength cover- 
age and S/N, we have designed the CORALS survey to have 
adequate spectral range and S/N to achieve a uniform lower 
redshift cut-off. 
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Fig. 3. The mass density of neutral gas, SIdla, in DLAs. 
Open circles and squares are measurements from the lat- 
est compilations by Peroux et al. (2001b) and Rao & 
Turnshek (2000) respectively. The solid circle is the value 
from the CORALS survey presented here for the redshift 
interval 1.8 < z a bs < 3.5. 



and 



b = - = 0.01, 

c 



(3) 



In cases where a Lyman limit system is present in the 
QSO spectrum within the redshift range considered, 



(zlls + 1) x 912 
1216 



(4) 



The total redshift interval for the CORALS sample 
is Az — 55.46 within which we detect 17 DLAs. Thus, 
n(z) = 0.31^5x08 at a mean absorption redshift (z Q b s ) = 
2.37 (the errors were calculated using the prescription by 
Gehrels 1986 for small number statistics). 

Taken at face value, the estimate of n(z) from the 
CORALS survey is ~ 50% larger (at the same {z a bs)) 
than that found in previous surveys. For example, Storrie- 
Lombardi & Wolfe (2000) deduced n(z) = 0.055(1 + 
z) 111 = 0.21 (no errors quoted). However, this difference 
is only marginally significant, since the two determinations 
of n(z) are within ~ la of each other. 



6.2. The Mass Density of Neutral Gas in DLAs 

The mass density of neutral gas in DLAs as a fraction of 
the closure density is expressed as 



DLA 



H Lim H 



CPcrit JN 



Nf(N)dN 



(5) 



where // is the mean molecular weight (= 1.3), run is 
the mass of a hydrogen atom, and f(N) is the column 
density distribution. We avoid making a priori assump- 
tions about the functional shape of f(N) by adopting the 
approximation by Storrie-Lombardi, McMahon & Irwin 
(1996) 



N„ 



Nf(N, z)dN 



AX 



9. 



so that 

H ol xm H Ei^Vi(H I) 



DLA 



CPcrit 



AX 



(6) 



(7) 



from which the mass density of neutral gas traced by 
DLAs is obtained by direct summation of the individual 
values of neutral hydrogen column density. The redshift 
path, AX, which takes into account co-moving distances, 
is given (in our adopted cosmology) by 



AX = £-[(! + *, 



i) a -l] (8) 



From eqs. (7) and (8) with the values of N(H I) listed 
in Table 3 we deduce logf^DLA^ = ~^-^-o.m (with errors 
calculated as described by Storrie-Lombardi et al. 1996). 
As can be seen from Figure ||, we again find that the 
CORALS value of Odla is in good agreement with previ- 
ous determinations. Thus, the most straight-forward con- 
clusion from the CORALS survey is that existing magni- 
tude limited samples of QSOs do not seriously underesti- 
mate the number of DLAs, nor their overall mass content; 
by and large, they seem to provide a fair census of neutral 
gas at high redshift. On the other hand, given the current 
statistical uncertainties, the data still admit a moderate 
degree of dust bias. The la limits on both n(z) and SIdla 
include the possibility that analyses based on optically se- 
lected samples may have underestimated both quantities 
by a factor of ~ 2 . This possibility is further highlighted 
when we consider the distribution of values of AT(H I) in 
the CORALS survey (Figure |j). Given the small size of 
our survey, we may well be missing DLAs at the high col- 
umn density end of the distribution simply through small 
number statistics. We return to this point in section 6.3 
below. We note that although it is impossible to constrain 
the H I column density distribution for the CORALS sam- 
ple, we do detect 2 relatively rare high N(H I) systems. 
Therefore, although we do not find evidence for a pre- 
viously excluded population of very high column density 
absorbers, with improved statistics it would be of great 
interest to re- assess the N(H I) distribution function. 
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Log N(HI) 

Fig. 4. Distribution of H I column densities in the large 
sample of Storrie-Lombardi & Wolfe (2000) and in the 
CORALS survey (shaded histogram). Due to the relatively 
small number of QSOs, and therefore DLAs, included in 
CORALS, we do not fully sample the iV(H I) distribution. 

We also consider another effect. As explained in section 
2, many previous surveys have included in their statisti- 
cal analyses candidate DLAs, identified on the basis of the 
equivalent width of the Lyman a line rather than by pro- 
file fitting to a damped profile. It is worthwhile examin- 
ing the overestimate of JIdla which may result from this 
approximation. For the 17 DLAs for which we have ob- 
tained our own spectra, we measure the equivalent width 
and compare the implied column density to that deter- 
mined by fitting the Lyman a line. In most cases we 
find the two techniques to be in very good agreement, 
certainly within the errors associated with each method. 
There are only three exceptions (the DLAs in B1055— 301, 
B1251-407a and B2314-409) where the equivalent width 
determination leads to a much higher A^H I) than the 
line fit. Inspection of Figure |^ shows that this is due to ex- 
tended absorption around the DLA. Although fits of these 
DLAs were not straight-forward, this process was facili- 
tated by higher spectral resolution and coverage of metal 
lines which provide additional guidance in the shape of 
the wings and z a b s - Had we used the values of N(H I) de- 
duced from the equivalent widths for the entire CORALS 
sample, we would have over-estimated ^dla by 20% . This 
discrepancy would have been further increased if extended 



blends of lines which do not include a DLA, such as those 
present in the spectrum of B1251— 407, were mistakenly 
included. Nevertheless, the over-estimate is not large and, 
given the increasing body of accurate measurements of 
A^H I) in DLAs, we think it very unlikely that this effect 
could be masking a higher degree of dust bias than that 
indicated by inspection of Figure 3. 

Finally, for completeness, we calculate ^dla in the 
redshift range 3.5 < z a bs < 4.0 where we detect two 
DLAs despite the fact that with CORALS we only sam- 
ple a total interval Az = 1.26 . The error bars are nat- 
urally very large, but all the same it is intriguing that 
logSlDLA^i seems to remain high at — 2.37L f Q' 2 g in con- 
trast with the slight down-turn suggested by the work 
of Storrie-Lombardi & Wolfe (2000) and Peroux et al. 
(2001b). It will be very interesting to see how better statis- 
tics will impact upon the CORALS value of SIdla at the 
highest redshifts since the present determination would 
suggest an increase in the importance of dust bias with 
increasing redshift. 

6.3. DLA Statistics as a Function of B -band 
Magnitude 

In their preliminary analysis of this sample, Ellison et al. 
(2000) found tentative evidence that ^dla was higher to- 
wards fainter QSOs, consistent with the effect expected 
from a dust bias. We re-examine this point in Figure ^, 
which shows cumulative statistics for the CORALS DLAs 
as a function of the i?-band magnitude of the background 
QSOs. Since there are relatively few bright QSOs in our 
sample, we also show the statistics for the LBQS DLA 
survey (Wolfe et al. 1995) which has a limit B ^ 19, using 
the column densities reported by Wolfe et al. (1995) and 
Storrie-Lombardi & Wolfe (2000). We confirm the initial 
conclusion by Ellison et al. (2000) that IIdla increases 
as fainter QSOs are observed, but stabilises at B ~ 20.0. 
Thus, we do not find a population of high N(H I) DLAs 
which is only revealed when faint (B ^ 20) QSOs are ob- 
served. However, closer inspection of the data emphasises 
the need to extend our survey in order to fully sample the 
column density distribution function, particularly at the 
high column density end. The increase in ^dla between 
the B < 19.0 and B < 19.5 bins is almost entirely due to 
a single DLA with N(H I)= 3.5 x 10 21 cm" 2 . Similarly, 
the increase between the B < 19.5 and B < 20.0 bins is 
caused by a single DLA with N(R I)= 4.5 x 10 21 cm~ 2 . 
Between them, these two systems account for over half of 
the neutral gas in the entire DLA sample. A K-S test that 
compares the distribution of N(H I) among DLAs towards 
B < 20 QSOs with those from the large sample of Storrie- 
Lombardi & Wolfe (2000) provides inconclusive results, 
i.e. that the samples are indistinguishable at only the la 
level. Taken as a whole, the CORALS DLAs are inconsis- 
tent (again only at the la level) with the N(H I) distri- 
bution of Storrie-Lombardi & Wolfe (2000). As found by 
Ellison et al. (2000) for a sub-sample of the CORALS sam- 
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Fig. 5. DLA statistics as a function of magnitude. The 
quantities plotted are for the cumulative values in the bins 
< the _B-band magnitude on the rc-axis. Bottom panel: 
Cumulative histogram of number of DLAs. The dashed 
area shows the number of LBQS DLAs used to determine 
improved statistics towards bright QSOs. Middle and top 
panels: stars show the values for CORALS DLAs, whilst 
the solid square shows the same quantity for the LBQS 
sample (B ^ 19); this point has been offset on the x-axis 
for clarity). 

pie, we confirm that for the complete sample that there 
is no strong correlation between magnitude and redshift 
out to z ~ 3.5. Therefore, these trends are not likely to 
be associated with the evolution (apparent or real) of the 
properties of the QSOs themselves. 

In addition to the cumulative values of n{z) shown 
in Figure ||, we calculate the number density of DLAs 
towards QSOs with B > 20 and B < 20, and find 
n(z) = 0.38t°; 2 ° (at (z abs ) = 2.54) and 0.27±g;J£ (at 
(zabs) = 2.25) respectively. For the B < 20 subset, this 
value is consistent with the number density found by 



Storrie-Lombardi & Wolfe (2000). Again, we see that there 
is an excess of DLAs in faint QSOs, but only at the ~ la 
significance level. 

7. Summary and Discussion 

We have presented the first results from the CORALS 
survey for DLAs in a radio-selected sample of QSOs. The 
sample consists of 66 z cm > 2.2 QSOs, 58 of which have 
been observed by us using the ESO 3.6 m, the AAT and 
the VLT facilities, while the remaining eight were culled 
from the literature. All the new spectra are presented in 
Figure 1. A total of 22 DLAs has been identified; 19 of 
which are at intervening redshifts. 

We find that the comoving mass density of neutral 
gas implied by these DLAs is logflDLA^ = — 2-59^ol24 
at a mean (z a bs) = 2.37, in good agreement with previ- 
ous surveys. Similarly, the number density of DLAs per 
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unit redshift in radio-selected QSOs, n(z) = 0.31 
only ~ lcr higher than that in optical, magnitude-limited 
samples. Within our own sample we also find that n(z) is 
higher, but again by only ~ lcr, in QSOs with B > 20, 
compared with sightlines towards brighter quasars. These 
results indicate that at redshifts z = 2 — 3.5 DLA surveys 
using optically selected QSOs probably underestimate the 
number of DLAs, and the gas mass they trace, by no 
more than a factor of about 2, in broad agreement with 
the predictions by Pei & Fall (1995). In particular, we 
have not uncovered a population of high column density 
A^H I)> 10 21 cm -2 absorbers which had been missed in 
previous searches limited to QSOs brighter than B ~ 20. 

These conclusions are somewhat tentative because of 
the small size of the CORALS sample. Our value of SIdla 
is dominated by two very high column density systems, 
both of which occur in moderately bright (B = 19.5,20) 
QSOs, and the column density distribution function of 
DLAs is clearly not well sampled with the relatively small 
number of QSOs in our survey. 

The next important step in this work is to determine 
the metallicities and dust content of CORALS DLAs. The 
possibility of a dust bias in DLA selection has previ- 
ously been appealed to in order to explain the lack of 
high column density, metal-rich absorbers (e.g. Prantzos 
& Boissier 1999). This explanation now seems less likely, 
regardless of whether or not CORALS DLAs prove to 
be more metal-rich, simply because high column density 
DLAs do not appear to be significantly more common in 
fainter QSOs, at least within the statistical limitations 
of our survey. Nonetheless, determining the metallicities 
of the new DLAs discovered here remains an important 
goal, because it will allow us to assess whether the low 
element abundances found so far are indeed typical of the 
full DLA population. 

What could then be the reason for the observed dearth 
of high column density, metal-rich DLAs? We consider it 
unlikely that gravitational lensing may be the answer. In 
principle one may conjecture that close alignment of QSOs 
with foreground galaxies may produce a tendency for such 
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sightlines to be deflected away from the inner regions of 
galaxies, where interstellar clouds with high N(H I) and 
high Z may be preferentially intercepted. However, to date 
no statistical evidence for lensing of QSOs by DLAs has 
been found (Le Brun et al. 2000) and, in any case, lens- 
ing would be most effective at significantly lower redshifts 
than those considered here (Smette, Claeskens & Surdej 
1997). 

A more plausible explanation is that there is simply a 
cross-sectional bias against detecting DLAs in sightlines 
that pass through the centres of galaxies. Observationally, 
one could argue that high N(K I), high Z DLAs have al- 
ready been found in the Lyman break galaxies (e.g. Pettini 
et al. 2000) which indeed have typical linear sizes one 
order of magnitude smaller than the impact parameters 
of most DLAs (Giavalisco, Steidel, & Macchetto 1996; 
Calzetti & Giavalisco 2000; see also Figure 1 of Pettini 
2001). Theoretical studies of DLAs also support this in- 
terpretation. For example, the models of Mathlin et al. 
(2001), who simulated DLA surveys by sampling model 
galaxies at random impact parameters, predict that the 
locus of high column density and metal rich absorbers 
should be populated, but that DLAs with these proper- 
ties are intrinsically rare due to the small cross-sectional 
area presented by the inner galactic regions where they 
are found. 

As already emphasised, the CORALS data set is too 
small to sample properly the column density distribution 
function and a considerably larger survey is required in or- 
der to provide the statistical coverage that will determine 
the true incidence of A^(H I)> 10 21 cm -2 absorbers. A 
statistically larger survey will not only improve our deter- 
minations of ^dla and n(z), but will also make it possible 
to investigate the evolution (or lack thereof) these quan- 
tities with redshift. It will be particularly interesting to 
examine the possibility of higher dust bias at larger red- 
shifts as suggested by the data in Figure 3; such data will 
offer a n insight into the evolution of dust at early epochs , 
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One prl omising prospect for extending the work presented 
here is the FIRST QSO survey (Gregg et al. 1996; White et 
al. 2000). Although spectral follow-up has so far been lim- 
ited to bright R < 19 targets, future follow-up of FIRST 
sources at fainter optical magnitudes would provide an 
excellent complement to the CORALS survey. 

Acknowledgements. It is a pleasure to acknowledge the con- 
sistent support of this project by the ESO and AAT Time 
Assignment Panels and the professional and efficient help of 
the telescope staff at the AAT, ESO 3.6-m and VLT. tn par- 
ticular, we are grateful to the ESO Paranal science operations 
staff for their expert execution of our service observations. We 
thank Mauro Giavaiisco and Lisa Storrie-Lombardi for obtain- 
ing spectra of two of our targets and Joop Schaye for usefuf 
comments on an earfier draft of this paper. 



ph/99112041 ) 



References 

Boisse, P., Le Brun, V.. 
1998, A&A, 333,841 



Bergeron, J., & Deharveng, J.-M., 



Pei, Y., & Fail, S. M., 1995, ApJ, 454, 69 

Pei, Y., Fall, S. M., & Bechtold, J., 1991, ApJ, 378, 6 

Peroux, C, Storrie-Lombardi, L. J., McMahon, R. G, Irwin, 

M., & Hook, I. M., 2001a, AJ, 121, 1799 
Peroux, C, McMahon, R. G, Storrie-Lombardi, L. J., & Irwin, 

M., 2001b, MNRAS, submitted 
Pettini, M. 2001, in The Promise of First, eds. GL. Pilbratt, 

J. Cernicharo, A.M. Heras, T. Prusti, & R. Harris, ESA 

SP-460, in press 
Pettini, M., Ellison, S. L., Steidel, C. C, & Bowen, D. V. 1999, 

ApJ, 510, 576 

Pettini, M., Smith, L.J., King, D.L., & Hunstead, RW. 1997, 
ApJ, 486, 665 

Pettini, M., Steidel, C.C., Adelberger, K.L., Dickinson, M., & 

& Giavalisco, M. 2000, ApJ, 528, 96 
Prantzos, N., & Boissier, S., 2000, MNRAS, 315, 82 
Prochaska, J. X., & Wolfe, A. M. 1999, ApJS, 121, 369 
Prochaska, J. X., & Wolfe, A. M. 2000, ApJ, 533, L5 
Rao, S.M., & Turnshek, D.A. 2000, ApJS, 130, 1 
Sargent, W., Steidel, C, & Boksenberg, A., 1989, ApJS, 69, 

703 



S. L. Ellison et al.: CORALS I: DLA number density and 



Savaglio, S., 2000, Invited talk at the IAU Symposium 204 'The 
Extragalactic Infrared Background and its Cosmological 
Implications', M anchester, August 2000, eds. M. Harwit 



& M.G. Hauser, a,stro-ph/0011473 



Shaver, P. A., Wall, J. V., Kellermann, K. I., & Jackson, C. 

A., Hawkins, M. R. S., 1996, Nature, 384, 439 
Songaila, A. 1998, AJ, 115, 2184 

Smette, A., Claeskens, J.-F., & Surdej, J., 1997, New 

Astronomy 2, 53 
Storrie-Lombardi, L., Irwin, M., & McMahon R. G., 1996, 

MNRAS, 282, 1330 
Storrie-Lombardi, L., McMahon R. G., & Irwin, M., 1996, 

MNRAS, 283, 79 
Storrie-Lombardi, L., & Wolfe, A. M., 2000, ApJ 543, 552 
Tytler, D., Fan, X.-M., & Buries, S. 1996, Nature, 381, 207 
Turnshek, D. A., Wolfe, A. M., Lanzetta, K. M., Briggs, F. H., 

Cohen, R. D., Foltz, C. B., Smith, H. E., & Wilkes, B. J., 

1989, ApJ, 344, 567 
White, R., et al, 2000, ApJS, 126, 133 

Wolfe, A. M., 1988, QSO absorption lines: Probing the uni- 
verse; Proceedings of the QSO Absorption Line Meeting, 
296 

Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 

1986, ApJS, 61, 249 
Wolfe, A. M., Lanzetta, K. M., Foltz, C. B., & Chaffee, F. H., 

1995, ApJ, 454, 698 



14 



S. L. Ellison et al.: CORALS I: DLA number density and $1dla- 



B0017-307 z = 2.666 



B0039-407 z = 2.478 



3x10" 



2x10" 



10" 




300 



200 



100 




4000 



5000 



3500 4000 4500 



B0104-275 z = 2.492 



B0113-283 z = 2.555 



5x10 



4x10 



2x10 




150 



100 




4000 



5000 



3500 4000 4500 



B0122-005 z = 2.280 

1000 H ' 



B0244-128 z = 2.201 




3500 4000 4500 
Wavelength 




3500 4000 4500 
Wavelength 



Fig. 1. Spectra obtained for the CORALS QSOs at the AAT (RGO spectrograph), ESO 3.6-m (EFOSC), and VLT 
(FORS1) facilities. Note that neither the higher resolution AAT spectra (FWHM = 3 A) nor the VLT ones are flux 
calibrated and the j/-axis is simply relative counts uncorrected for instrument response and atmospheric transmission. 
The ESO 3.6-m data (FWHM ~ 7 A), however, were flux calibrated and the y-axis scale is in units of ergs s 1 cm 2 
A -1 . Vertical tick marks show the position and redshift of intervening DLAs, if present. 
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Fig. 1. Continued 
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Fig. 1. Continued 
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Fig. 1. Continued 
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Fig. 1. Continued 



20 



S. L. Ellison et al.: CORALS I: DLA number density and J7dla- 



B1354-107 z = 3.006 



B1402-012 z = 2.518 




8x10" 



6x10" 



4x10" 



2x10" 




3500 4000 4500 



4000 



5000 



B1406-267 z = 2.430 



B1418-064 z = 3.689 



600 



400 



200 





4000 



5000 



5000 5500 6000 



B1430-178 z = 2.331 



B1535 + 004 z = 3.497 



4x10 



-16 



2x10" 




4000 5000 
Wavelength 



400 



200 




4000 5000 
Wavelength 
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Fig. 1. Continued 
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Fig. 2. Damped Lyman a profile fits (continuous lines) to all newly discovered intervening DLAs. See Table 3 for the 
values of N(H I) and z a b s corresponding to the theoretical profiles shown. Not that the bottom right-hand panel is on 
a different wavelength scale due to the large column density of the DLA towards B1055— 301. 
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